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Sh ' Abstract. A class of solvable (systems of) nonlinear evolution PDEs in multidimensional 

space is discussed. We focus on a rotation-invariant system of PDEs of Schrodinger type and 
on a relativistically-invariant system of PDEs of Klein-Gordon type. Isochronous variants 
of these evolution PDEs are also considered. 
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1 Introduction 



Over a decade ago a class of C-integrable - i.e., solvable via a Change of variables - systems 
of PDEs in multidimensional space were identified [2]. (A problem involving nonlinear PDEs 
is considered solvable if its solution can be obtained by performing algebraic operations - such 
as finding the zeros of a given polynomial - and by solving linear PDEs; of course only seldom 
these operations can be performed explicitly.) In the present paper - motivated by the scarcity 
of solvable models of nonlinear evolution PDEs in multidimensions hence by the interest of any 
such model - we study (a subclass of) these solvable PDEs in more detail than it was done 
hitherto. We focus mainly on the system of PDEs of Schrodinger type 

•/ a/ , urr-\ i o a + byj n yj m - (Vtp n ) ■ (Vyj m ) 

m=l, mj^n 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue "Integrable Systems and Related Topics" 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 
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which is rotation-invariant if W(f) = W(r), and on the relativistically-invariant system of PDEs 
of Klein-Gordon type 



N 

a 



V'n.tt - AV>„ + M = 2 2^ / _ / • I 1 " 2 ) 

m=l, m^=n 

Notation. Here and throughout N is an arbitrary positive integer (N > 2); the index n, 
as well as other analogous indices (see below), range generally from 1 to N; the dependent 
variables ip n = ip n (r, t) are generally considered complex (although this is only mandatory for 
the first, (jl.lh . of these two systems of PDEs); the space variable r is a (real) ^-vector (with 
S arbitrary: for S = 3, r = (x,y,z)), and r indicates its modulus, r 2 = r • r; a dot sandwiched 
among two S'-vectors denotes the standard (rotation-invariant) scalar product (for instance, 
for S = 3, r-f=x 2 +y 2 + z 2 ); the "potential" W (r) is an arbitrary function of the spatial 
coordinate r (we will often assume that it only depends on the modulus r of the S- vector r); the 
constants a and b, as well as the "mass" parameter M, are arbitrary (they might also vanish); t is 
the (real) time variable; subscripted independent variables always denote partial differentiations; 
V respectively A = V • V are the gradient respectively the Laplace operator in 5-dimensional 
space. Isochronous variants of these evolution PDEs are also considered. 

In the following Section 2 we review tersely the general class of PDEs treatable in this 
manner and the technique to solve them. In Sections 3 respectively 4 we treat in some detail 
the PDEs (jl.ip respectively (jl.2p . describing various properties of their solutions and reporting 
some representative examples, and we exhibit their isochronous versions. In Section 5 we take 
advantage of the electronic format of this article to present a few animations displaying visually 
the time evolution of a few solutions of some of these solvable (systems of) nonlinear evolution 
PDEs: in this article we restrict these presentations to very few cases, all with space dimensiona- 
lity less than three (we are of course aware that the three-dimensional case is probably the most 
interesting one - since we seem to live in a three-dimensional world - but the presentation of 
animations in a three-dimensional context is somewhat more tricky and we therefore postpone 
the display of such examples to future articles we hope to issue soon). In Section 6 we outline 
future directions of research. 



2 A class of solvable (systems of) PDEs 
in multidimensional space 

In this section we review tersely the basic idea allowing to identify a class of solvable nonlinear 
evolution PDEs. The interested reader will find a more detailed treatment in the paper where 
this approach was introduced [T], and especially in [3] where this method is treated in con- 
siderable detail: see Section 2.3 of this book, and other references quoted there in Section 2.N. 
These treatments focussed however on ODEs rather than PDEs: the extension to PDEs is rather 
straightforward, although it took some time to realize its feasibility [2] (see also Exercise 2.3.4.2-5 
in $)■ 

Let ^ (ip\ r, t) be a (time- and space-dependent) monic polynomial of degree N in the vari- 
able tp, and denote by ip n (f, t) its N zeros and by ip m (f, t) its N coefficients: 

N 

m(^?t) = J[[^-i> n (r,t)], (2.1a) 

n=l 

N 

* (V; r, t) = ^ + <Pm (r, t) ^ N ~ m . (2.1b) 

m=l 
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Assume then that the time evolution of the dependent variable ^ - hence its dependence on 
the time and space independent variables t and f, as well as its dependence on the independent 
variable ip - is characterized by a linear evolution PDE, which must of course be consistent with 
the fact that ^ is, for all time, a monic polynomial of degree N in ifj, see (|2.ip . The linear 
character of this PDE entails that the time-evolution of the N coefficients ip m (r, t), see (|2. lb|) . 
is as well linear, possibly explicitly solvable (see below). On the other hand the corresponding 
time evolution of the N zeros ip n (r, t) will be nonlinear, due to the nonlinear character of the 
mapping relating the iV zeros tp n (r, t) to the N coefficients ip m (r, t) of the polynomial (tp; f, t), 
see (|2.ip . These (systems of) nonlinear evolution PDEs satisfied by the dependent variables 
V'n ifi t) are those referred to in the title of this paper: they are indeed generally solvable by 
taking advantage of the mapping, see (12. II) . relating the N dependent variables ip n (r,t) to the 
N functions ip m (r , t) . 

In particular it is easily seen (using, if need be, the formulas provided in Section 2.3 of [3]), 
that to the evolution PDE 

W t _ A* - V(r) (V^v " NV) + o*^ + b [iP 2 *^ - N (N - 1) *] = (2.2) 

satisfied by ^ (ijj;r,t) , there corresponds for the N zeros ijj n (f,t) , see (|2.1aj) . just the system 
of nonlinear evolution PDEs of Schrodinger type with 

W (f) = V(r) - 2 (N - 1) b, (2.3) 

while the corresponding evolution of the N coefficients (p m (r,t), see (|2.1bp . is clearly given by 
the system of linear evolution PDEs 

iipm t t - Aip m + [W (f) - b (m + 3)] mip m 

= -a(N -m + 2)(N -m + l)^ m _ 2 , m = l,...,N (2.4) 

with if-\ = and ifo = 1 (see (I2.1bj) ). Note that this system of linear PDEs is decoupled if the 
constant o vanishes; as we shall see in the following section, it can also be replaced by a decoupled 
system if the constant a does not vanish, but the potential W (r) is constant, W(r) = C, 

see Section 3. 

Likewise to the evolution PDE 

tt« -Af-/i 2 [V^ " NV] + o*^ + b [V 2 ^^ - N ( N ~ 1) *] = °» ( 2 - 5 ) 

there corresponds for the iV zeros ip n (r, t) just the system of nonlinear evolution PDEs of 
Klein-Gordon type (II. 2p with 

M 2 = fj 2 — 2(N — 1) b, 

while the corresponding evolution of the iV coefficients (p m (r, t) is clearly given by the system 
of linear evolution PDEs 

Aip rn + [M 2 -b(m + 3)] ?ntp m 
= -a{N -m + 2){N -m + l)(p m _ 2 , m = l,...,N, (2.6) 

again with cp_i = and tfo = 1 (see (j2.1b|) ): this system is decoupled if the constant a vanishes, 
and can be replaced by a decoupled system even if a does not vanish, see Section 4. 

Having conveyed, tersely but hopefully clearly, the main idea of this approach to identify 
solvable systems of nonlinear evolution PDEs, we turn, in the next two sections, to the study 
of the two systems of nonlinear evolution PDEs (11. ip and (II. 2p . 
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3 Solvable system of nonlinear PDEs of Schrodinger type 

In this section we investigate the system of nonlinear evolution equations (jl.lj) , firstly by analytic 
techniques and subsequently by reporting some of its solutions in a representative set of cases. 

Remark 3.1. The "mean field" tp{r,t) 
1 N 

^) = ^I>n(r,t) (3.1) 

n=l 

satisfies the linear Schrodinger equation 

i%pt - Aip + W (r) tp = 0. 

Proof. Sum the nonlinear evolution PDEs (jl.ip over n from 1 to N, use (|3. 1 j) in the left-hand 
side, and notice that the double sum in the right-hand side vanishes due to the antisymmetry 
of the summand under the exchange of the two dummy indices n and m. ■ 

Remark 3.2. If all the N coefficients tp m vanish, ip m = 0, then = ip N (see (|2.1b|) ). hence all 
its N zeros ip n correspondingly vanish, tp m = (and more generally: if the first M coefficients (p m 
vanish, (p m = for m = 1, . . . , M, then M of the N zeros of \I/ vanish). Hence to a set of 
localized solutions ip n (r , t) of the system of nonlinear evolution equations (jl.ip , characterized 
by the asymptotic conditions 

lim [Vn (r, t)]=0, n = 1, . . . , N, (3.2a) 



there correspond a set of localized solutions of the system of linear Schrodinger PDEs (12. 4h 
(with (|2,3p ). characterized by the analogous asymptotic conditions 

lim [tpn(r,t)} =0, n = l,...,N; (3.2b) 

1 — >oo 

and, of course, viceversa, namely clearly (|3.2b|) entails (|3.2ap . However when two different 
components, ip n (r, t) and tp m (f, t) with n / m, of the system of nonlinear evolution equa- 
tions (jl.ip are equal, ijj n (r,t) = ip m (r,t), clearly this system runs into a singularity due to the 
vanishing of one of the denominators in the right-hand side of (jl.l|) ; and this situation gets 
even worst if the two different components, ip n (r, t) and ip m (r, t) with n ^ m, both vanish, 
ip n (r, t) = ipm (f, t) = 0. Hence in all the examples considered below we shall try and avoid this 
problem, and in particular we shall focus, rather than on localized solutions ip n (r, t) that vanish 
asymptotically (see (|3.2ap ). either on solutions that oscillate asymptotically, or on solutions that 
are asymptotically constant, 

lim hjj n (r,t)] = a n with a n ^ a m if n ^ m. (3-3) 

r— >oo 

Note that these asymptotic values a n might depend on the direction along which the space 
coordinate r diverges. 

Remark 3.3. If the potential V (r) is constant, 

V(r) = B, (3.4) 
(entailing that W (r) is as well constant, 



W (r) = B - 2 (N - 1) 6, 
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see (|2.3p ). it may be convenient to replace the expression (|2.1bp of the monic polynomial 
^ (ift; f, t) in terms of its coefficients ip m (r, t) by the following representation: 

* ^ r ' t} = Ib^T + £ Xm (r ' } ( a) 

fe\ 1/2 

--J , (3.5b) 

where the polynomial C] (z), of degree £, is the standard Gegenbauer polynomial [5], satisfying 
the ODE 

(1 - z 2 ) (z) - (2 7 + 1) zC?" (z) + £(£ + 2 7 ) (*) = 

(where appended primes denote derivatives with respect to the argument of the function they 
are appended to, in this case with respect to z), and being characterized by the asymptotic 
behavior 



lim 



& (z) 



2*r( 7 + i) 

1, where k e = ^ ^ . 



Then the requirement that \l/(^>;r, i) satisfy the PDE (|2.2p with (|3.4p . and therefore that its 
TV" zeros ip n (f, t) satisfy the system of nonlinear PDEs with (|3.4p , entails that the iV 

coefficients Xm (r, t) satisfy the (system of decoupled) linear Schrbdinger PDEs (with constant 
potentials) 

i Xm ,t ~ A Xm + [B - b(2N -m- 1)] m Xm = 0, m = 1, . . . , N. (3.6) 

The investigation of the limiting cases in which some coefficients vanish, and therefore a different 
set of classical polynomials come into play in place of the Gegenbauer polynomials, is left as an 
exercise for the diligent reader. 

The nonlinear mapping among the N dependent variables ip n (r, t) satisfying the system of 
nonlinear PDEs (jl.ip and the N functions <p m (r, t) respectively %m (^ t) satisfying the system of 
linear PDEs (|2.4p respectively (|3.6p . entailed by the simultaneous validity of (|2.1aj) and (|2.1bj) 
respectively (|3.5|) . is the key to the solvability of the system of nonlinear PDEs (jl.ip . This 
can be taken advantage of in two ways: to solve the initial-value problem for the system of 
nonlinear PDEs (jl.ip . or to manufacture special, possibly quite explicit, solutions of this system 
of nonlinear PDEs. 



3.1 Solution of the initial-value problem for the system 
of nonlinear Schrodinger-like PDEs ( 11. lj) 

The initial-value problem consists in the determination of the solution ifj n (r,t), n = 1,...,N, 
of the system of Schrodinger-like nonlinear PDEs (jl.ip corresponding to assigned initial data 
ip m (r,0),m = l,...,N. 

The first step is to determine the corresponding initial data (p m (r, 0) of the system of 
Schrodinger-like linear PDEs (j2.4p . This is achieved by solving for the N functions tp m (r, 0) the 
system 

N N 

i; N +J2Vm (f, 0) ^ N ~ m = J] bl> ~ V>n (F, 0)] , (3.7) 
m=l n=l 
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entailed by the simultaneous validity of (|2.1ap and (j2.1bp at t = 0. This amounts to the 
determination of the N coefficients ip m of a monic polynomial given its N zeros ip n ; the relevant, 
explicit formulas are of course well-known: 

N N 

V9i(f,0) = -^«(r,0), <p2(r,0)= £ ^ (f, 0) ^ m (r, 0) , 

n=l ti,m=l,m^n 

and so on, up to 

N 

<PN (r,0) = (-) 7V n^ 

n=l 

The second step is to solve the initial- value problem for the system of evolution PDEs (|2.4p . 
obtaining thereby its solution ip m (r, t) at time t. The linear character of this (coupled) system 
of evolution PDEs, (12. 4ft . provides the main simplification; of course an explicit solution is only 
possible for special choices of the potential W(r). 

The third step is to obtain the solution ifj n (r, t) from the, now assumedly known, functions 
ip m (r , t) , via the relation 

N N 

IJ - ^ (r, t)] =i, N +Y,Vm (r, t) i, N ~ m , (3.8) 
n m=l 

again entailed by the simultaneous validity of (I2.1al) and (|2.1b|) . but now at time t. This amounts 
of course just to the purely algebraic task of finding the zeros of a given monic polynomial of 
degree N: an explicit solution is of course only possible for iV < 4. 

In the special case of a constant potential V (r), see (|3.4p . an alternative procedure of solution 
can be based on the representation (|3.5p rather than ()2.1b[) : this eases the second of the steps 
outlined above, but makes a bit less simple the first step. The diligent reader will easily figure 
out the relevant details. 

3.2 How to manufacture explicit solutions of the system 
of nonlinear Schrodinger-like PDEs ( 11. lj) 

Clearly the appropriate strategy - underlining all the examples exhibited below - is to identify 
an explicit solution (p m (f,t) of the system of linear evolution PDEs (|2.4p . and then to obtain 
the corresponding solution ip n (f,t) of the system of nonlinear evolution PDEs (jl.lj) via (|3.8|) . 
namely by identifying the N zeros ip n (r, t) of the monic polynomial of degree TV in ip having 
the coefficients ip m (r,t), see (l2.1bp . This can of course be done explicitly only for N < 4: not 
a significant restriction when it comes to the explicit exhibition of examples, which would indeed 
be impractical for larger values of N (in Section 5 we indeed limit our exhibition of animations 
to the 3-body case, N = 3). 

An alternative route - applicable when the potential is constant, see ()3.4p - takes as starting 
point an explicit solution Xm if it) of the [decoupled) system of linear evolution PDEs (|3.6p . and 
then obtains the corresponding solution tp n (r, t) of the system of nonlinear evolution PDEs (11. ip 
via ()3.5p . namely by identifying the N zeros (r, t) of the monic polynomial of degree N in ijj 
given by this expression (|3.5p . 

3.3 Examples 

Example 3.3.1. The simplest example is characterized by the assignment 



a = b = 0, W(r) = V(r) = 0, 
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namely by the system of nonlinear evolution PDEs in one-dimensional space (see (jl.ip ) 

Jlpn ■ V' 



■I A I O ^« ' ^ m /o ^^ 

^n,t-A^ n = -2 ^ _,/, ■ ( 3 - 9 ) 



m=l, rn=£n 

The corresponding system of linear PDEs satisfied by the coefficients (p m (r, t) of the poly- 
nomial (ip] f, t) in the variable ip of which the solutions ip n (f, i) are the iV zeros, 

N N 

* (V; r, t) = ^ + 2 *) ^ 7V_m = II & ~ ^ *)] • ( 3 - 10 ) 

m=l n 



reads as follows (see (|2.4p ): 

i(p m>t - Aip m = 0, m = l,...,iV. (3.11) 

A special class of "traveling wave" solutions of these (decoupled) linear PDEs reads 

(r,t) =cp m (f- vt) = A m + B m exp [-iv(r- vt)] , (3.12) 

with v an arbitrary real constant 5- vector and the 2N scalar constants A m , B m also arbitrary 
(possibly complex). These solutions (which are clearly the most general ones of "traveling 
wave" character, namely depending on the single S- vector f — vt rather than separately on the 
S'-vector space variable r and the scalar time variable t) are not localized: they are constant 
along the S— 1 space directions orthogonal to v, periodic with period L = \2ir/v\ along the space 
direction parallel to v, and periodic in t with period T = 2ttv~ 2 . The corresponding traveling 
wave solutions tp n (x,t) of the system of nonlinear PDEs (|3.9I) are the N zeros of the following 
polynomial of degree N in the variable ip: 

N 

Pn WO + Qn-i (V>) exp [-iv(f- vt)] = JJ [ij, - V>„ (r,t)} , (3.13a) 

n 

where Pj\r (ip) is an arbitrary monic polynomial of degree N and Qjv-i (ip) is an arbitrary 
polynomial of degree N — 1, 

N N 

P N (V) = ^ + A m^ N ~ m , QN-1 W = E ^^~ m ' ( 3 - 13b ) 

m=l m=l 

This entails of course that these solutions are as well of the same traveling wave type, ip n (r, t) = 
ipn if ' ~ vt), with the same periodicity properties described above for (p m ift) = ip m (r — vt) - 
or possibly with time periods which are integer multiples of that of the coefficients ip m (r, t) due 
to the possibility that through the time evolution the zeros of the polynomial (|3.13bl) exchange 
their roles; for each component ip n (f, t) this integer multiple cannot of course exceed N, while for 
the entire solution {ijj n (f,t) ;n = 1, . . . , N} the period can be somewhat larger - but generally 
not too much [6]. 

A solution of this kind (with S = 2 and N = 3, and a specific choice of the remaining free 
parameters) is displayed as an animation in Section 5. 

Another special set of solutions of (13, lip - written below, for simplicity, for the one-dimensional 
case (5=1)- reads 



(p m (x, t) = A m H — — exp 

{t + t m + ir/ r 



lm) 



i (x x m -\- iTj m v m ) 
4 (t + t m + ir] 

in / 



(3.14) 
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where the constants A m and B m are arbitrary (possibly complex; but see below for some con- 
ditions on the constants A m ), the constants x m , t m and v m are also arbitrary but real, and the 
constants r] m are also arbitrary but positive, r} m > 0: it is indeed easily seen that these conditions 
are sufficient to guarantee that (p m (x, t) is nonsingular for all real values of the independent 
variables x and t and tend asymptotically to the constants A m : 

lim [ip m (x,t)] = A m . 

x— »±oo 

The corresponding solutions ip n (x,t) of the system of nonlinear PDEs (|3.9p are the N zeros of 
the following polynomial of degree N in the variable ip: 



N 



Br, 



(t + t m + ir] m ) 1/2 



cxp 



m=l 

N 



i (x X m ~\~ irj-m^rn) 
^{t + t m + irj 

m ) 



'0 



N-m 



(3.15) 



ri=i 



Note that this entails that these solutions satisfy the asymptotic property (|3.3|) . with the con- 
stants a n being the N zeros of the polynomial, of degree TV" in the variable a, 



N 



N 



+ 2> 

m=l 



N-m 



J (a - o„) . 



n=l 



Of course the constants ^4 m should be assigned so that these N zeros o n are all different, see (|3,3p , 
A solution of this kind (with N = 3, and a specific choice of all the remaining free parameters) 
is displayed as an animation in Section 5. 



3.4 Isochronous version of the class of nonlinear PDEs of Schrodinger type 

In this subsection we report, without much commentary, an "cj-modified" version of (a subclass 
of) the system of nonlinear PDEs of Schrodinger type (jl.lj) . which is characterized by the 
property to possess an ample class of solutions completely periodic in time with period 



T 



2 k 



(3.16) 



(the reason why this is so will be rather obvious from what follows; for more details on, and 
other examples of, isochronous PDEs obtained in an analogous manner see [7] and Chapter 6 

of a). 

We start from the subcase of the system of nonlinear PDEs of Schrodinger type (Jl.lj) with 
b = 0, W(f) = 0, 
and we set 



i>n if, t) 



iXut 



ip n (p,r), p = p(t) = e 2 r, t = t (t) 



IUJ 



with lu an arbitrary positive constant and A an arbitrary real rational number if a vanishes 
(a = 0, see (jl.ip and below), otherwise A = — 1/2. It is then easily seen that the new dependent 
variables ip n (r, t) satisfy the following system of nonlinear PDEs of Schrodinger type: 



N 

•■ E 

m=l, mj^n 
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Note that this system is autonomous with respect to the time variable, but it features an explicit 
dependence on the space variable r (see the last term in the left-hand side); and it is clearly 
rotation invariant. 

Example 3.4.1. The simplest example is again characterized by the assignment 
5 = 1, O = o, 

namely by the system of nonlinear evolution PDEs in one-dimensional space (see (jl.ip ) 

i n.x 



n,x 



N 



m=l, m^n 

A special class of solutions of this system of nonlinear PDEs obtains, via the formula 

4> n (x, t) = exp (iXut) ip n (£, t) , 
£ = = exp 



iuit\ 
2 J 



r(t) 



exp (iujt) — 1 



(3.17) 



from the solutions ip n (£,t) of (|3.13ap . now reading 



N 



Pn W + Qn-i (VO exp [-iv (£ - vt)\ = JJ - ip n (?, r)] , 

n 

where Pn (tp) is an arbitrary monic polynomial of degree N and Qn-i (ip) is an arbitrary 
polynomial of degree N — 1, see (|3.13b|) . 

Another class of solutions is provided via (|3.17p from the solutions tp n of (|3.15p , now reading 



N 



m=l 

N 



(r + T m + ir]„ 
\[i>-^n (£,T)]. 



,1/2 



exp 



4 (r + T m + irj 



4' 



N-m 



,i = l 



4 Solvable systems of nonlinear PDEs of Klein-Gordon type 

In this section we investigate, firstly by analytic techniques and subsequently via the explicit 
display of a few of its solutions, the system of nonlinear evolution equations of Klein-Gordon 
type (|1.2p . The first remarks are analogous to those given in the first part of the preceding 
section and are therefore reported below without much commentary (their proofs are analogous 
to those given in the preceding section; we also use occasionally the same notation, confiding 
that this will cause no misunderstandings). 

Remark 4.1. The "mean field" tp (f,t) defined by (|3.ip satisfies now the linear Klein-Gordon 
equation 

4>tt - + M 2 4> = 0. 

Remark 4.2. To a set of localized solutions ip n (r, t) of the system of nonlinear evolution 
equations of Klein-Gordon type (11. 2h characterized by the asymptotic conditions (|3.2ap there 
correspond a set of localized solutions of the system of linear Klein-Gordon PDEs (|2.6p char- 
acterized by the analogous asymptotic conditions (j3.2b|) ; and, of course, viceversa. But such 
localized solutions cause the same kind of problem discussed in Remark 13.11 Hence in all the 
examples considered below we shall try and avoid this problem, just as indicated in Remark 13.21 
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Remark 4.3. It may be convenient to replace the expression (|2.1b|) of the monic polynomial 
^ (tp;f,t) in terms of its coefficients ip m (r,t) by the representation (|3.5a|) . again with (l375bl) 
but now with 

M 2 + b(2N-m) 
7 = 26 ' 

instead of (j3.5c|) . Then the requirement that (tp; r, t) satisfy the PDE (|2.5p . and therefore that 
its N zeros t/j n (r,t) satisfy the system of nonlinear PDEs (|1.2|) . entails that the N coefficients 
Xm (?■> t) satisfy the (system of) decoupled linear Klein-Gordon PDEs 

Xm.it " Ax™ + [M 2 - 6(m + 3)] m Xm = 0, m = 1, . . . , N. (4.2) 

Analogous developments to those reported in the preceding Section 3 in the context of the 
system of Schrodinger-like nonlinear evolution PDEs (jl.ip can now be elaborated in the present 
context of the system of Klein-Gordon-like nonlinear evolution PDEs (j 1 . 2|) . Our presentation 
below is more terse than in the preceding section, to avoid repetitions. 



4.1 Solution of the initial- value problem for the system 
of nonlinear Klein— Gordon-like PDEs (11.2ft 

The initial-value problem consists now in the determination of the solution ip n (r , t) of the system 
of Klein-Gordon-like nonlinear PDEs (|1.2p corresponding to assigned initial data ip m (r, 0) and 
ipm,t (r,0) , m = l,...,N. 

The first step is to determine the corresponding initial data <p m (r, 0) and <p m t (^ 0) of the 
system of Klein-Gordon-like linear PDEs (|2.6|) . This is achieved by firstly solving, as above, for 
the N functions (p m (r, 0) the system (|3.7p . and then by solving for the N functions (p m ,t (^0) 
the system 

N N ( N 

£ <p m>t (r, 0) ^ N ~ m = - E i (r, 0) [] W ~ ^ °)] \ ■ 

m=l n=l ^ m=l,m^n J 

The second step is to obtain the solution (p m (r, t) at time t of the system of evolution 
PDEs (|2.6p . Again, the linear character of this system of evolution PDEs provides the main 
simplification. 

The third step is, as above, to obtain the solution ip n (r, t) from the, now assumedly known, 
functions ip m (r,t) via the relation (|3.8p . amounting again just to the purely algebraic task of 
finding the N zeros of a given monic polynomial of degree N. 

An alternative procedure of solution can be based on the representation (|3.5ap (with (|3.5b|) 
and (14. lj) rather than (|3.5cp ): this eases the second of the steps outlined above, but makes a bit 
less simple the first step. The diligent reader will easily figure out the relevant details. 



4.2 How to manufacture explicit solutions of the system 
of nonlinear Klein— Gordon-like PDEs (11.2ft 

As above the strategy - that underlies all the examples discussed below - is to identify an 
explicit solution <p m (r,t) of the system of linear evolution PDEs (|2.6p . and then to obtain the 
corresponding solution ifj n (r, t) of the system of nonlinear evolution PDEs (jl.2p via (13. 8|) , namely 
by identifying the N zeros ip n (r, t) of the monic polynomial of degree iV in i/j with coefficients 
(p m (f,t), see (]2.1bp . 

An alternative route takes as starting point an explicit solution Xm (r,t) of the (decoupled) 
system of linear evolution PDEs (|4.2p . and then obtains the corresponding solution ip n (r,t) 
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of the system of nonlinear evolution PDEs (jl.2p by identifying the quantities ip n (r, t) as the 
N zeros of the monic polynomial of degree N in ip given by the expression (|3.5a|) with (|3.5b[) 
and (|iTTI) . 

4.3 Examples 

Example 4.3.1. The simplest example is characterized by the assignment 

a = 6 = 0, M = 0, 
namely by the system of nonlinear evolution PDEs (see (jl.2p ) 

N 

1pn,tVr 



^n,tt-^n = 2 ^ — . (4.3) 



m=l, m^n 

The corresponding system of linear PDEs satisfied by the coefficients ip m (r, t) of the poly- 
nomial ^ (ip; r, t) in the variable tp of which the solutions ip n (r, t) are the N zeros, see (|3,10p . 
reads as follows (see (|2.6p ): 

(p m>u - Ac^ m = 0, m = l,...,N. 

The general solution of these linear PDEs reads 

K 

<Pm(r,t) = J ^2f m k(r- u k t) , (4.4) 

k=l 

with the KN functions f m \~ (f) arbitrary and the K constant 5-vectors u& having unit length, 
Uk = 1, but being otherwise arbitrary. Of course these solutions ip m (r, t) are localized if the arbi- 
trary functions f m k (r) are themselves localized, but (motivated by Remark I4.2p we shall rather 
consider solutions that tend asymptotically to nonvanishing asymptotic values; and (p m (r, t) 
has the character of a traveling wave if K = 1. 

The corresponding solutions ip n (r, t) of the system of nonlinear PDEs (|4.3p are the N zeros 
of the following polynomial of degree N in ip: 

N K N 

* ty; r, t) = ^ + E ^ 7V ~ m E f"* ( f ~ = II ^ - ^ ( f > *)1 • 

m=l fc=l n 

Two solutions of this kind (with 5 = 2 respectively 5 = 1, N = 3, and a specific choice of 
the remaining free parameters) are displayed as animations in Section 5. 

Example 4.3.2. An analogous example - but reported here for simplicity in the two-dimensional 
case (5 = 2)- is characterized by the analogous assignment 

5 = 2, a = b = 0, M = 0, 

namely by the system of nonlinear evolution PDEs in two-dimensional space (see (jl.ip ) 

, _, _ 9 ^ Mm,t ~ Mm* ~ ^ m ,y 

WnM Yn,xx Wn,yy — " / . i i 

r^ 1 , Vn-Wm 
m=l, m^pn 

The corresponding system of (decoupled) PDEs satisfied by the coefficients <p m (x, y, t) reads 

(Pm,tt (Pm,xx (Pm,yy — 0' 
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A class of regular solutions of this system of PDEs reads 

(p m (x, y, t) = J (V (x - x m ) 2 + (y - y m ) 2 ^ [An cos(t) + B m sin(i)] + C m , (4.5) 

where A m , B m and C m are 3N arbitrary constants (possibly complex), and Jo(r) is the zeroth- 
order Bessel function of the first kind. A solution of this kind (with N = 3, and a specific choice 
of the remaining free parameters) is displayed as an animation in Section 5. 

4.4 Isochronous version of the class of nonlinear PDEs of Klein— Gordon type 

In this subsection we report, with even less commentary than in the (analogous) Subsection 3.4, 
an "^-modified" version of (a subclass of) the system of nonlinear PDEs of Klein-Gordon 
type (|1.2j) . which is again characterized by the property to possess an ample class of solutions 
completely periodic in time with period T, see f|3. 16|) . 

Now we start from the subcase of the system of nonlinear PDEs of Klein-Gordon type (11. 2h 
with 



6 = 0, 
and we set 

V>n (r, t) 



M = 0, 



ipn(p,r), P = p(t) 



T = r{t) 



ijj 



with uj an arbitrary positive constant and A an arbitrary real rational number if a vanishes 
(a = 0, see (|1 .2|) and below), otherwise A = —1. It is then easily seen that the new dependent 
variables ip n (r, t) satisfy the following system of nonlinear PDEs of Klein-Gordon type: 

ijj nt tt - A^ n - 2zw(r • V)i> ni t - i (2A + 1) uj^t 

- A (A + 1) uj 2 ij n - (2A + 1) oo 2 {r ■ V)$ n - uo 2 {r- V) 2 ^ 



N 

^ E 

m=l, m^n 



+ 



tpn ~ 1pm 



m,t 



i\ujijj m — iuj(r ■ V)if) n 



Note that this system is autonomous with respect to the time variable, but it features an 
explicit dependence on the space variable r (incidentally, in the last term in the left-hand side, 

s 

(r- V) 2 = (f • V)(r- V) = S(r-V)+ Yl r fc r iVfcVj, where of course rj denotes the j'-th component 

k,j=i 

of the S- vector r); and it is clearly rotation invariant. 



5 Animations 

In this section we show a few solutions of the Schrodinger and Klein-Gordon type problems 
treated in the previous sections, displayed as animations over time. 

Let us begin with a brief description of the methodology - implemented via a software written 
using Mathematica - employed to obtain the numerical results presented below. After assigning 
the space dimension S, the coefficients ip n (f,t) (see (|2.4p and (|2.6|) ) and the monic polynomial 
9(tp;r, t) of degree N in tp (see (|2.ip ). we compute the N zeros ip n (f,t) of this polynomial 
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Figure 5.1.1. The ip m (f,t) functions and the absolute values of the ip m (r,t) solutions for < t < 1/8. 

as follows: firstly we create in the (S + l)-dimensional space of the independent variables r, t 
a lattice; next, we use a root-finding routine to calculate, at an appropriately chosen point of 
the lattice, the (generally complex) values of the N zeros ip n (r, t) of the polynomial (ip; r, t), 
taking them in a generic order; and then we use an iterative root-finding procedure to calculate 
the N zeros of the polynomial \F (?/>; r, t) at any new point of the lattice, making use of the 
zeros previously calculated at the nearest points of the lattice so as to preserve the same initial 
ordering of the zeros. 

In the following we consider only examples with iV = 3. The animations are organized as 
arrays of synchronized subanimations, showing above the time evolution of the three ip m (r, t) 
functions, and below the corresponding time evolution of the three ip n (r, t) solutions. In each 
frame of the animation (namely, for a fixed value of the time variable), we display the values 
of the functions ip m (r,t) and ip n (r,t) (or their absolute values \ip m (r, t)\ and \ip n (r, t)\, when 
these functions are complex), with respect to the x variable in two-dimensional plots if S = 1, 
or with respect to the x and y variables in three-dimensional plots if S = 2. In this paper we 
restrict attention only to cases with space dimension S = 1 and S = 2. 

5.1 Solutions of ( TTTTI) 

In this subsection we present two numerical solutions of (jl.ip displayed as animations, the first 
one with S = 2 and the other two with S = 1. 

Example 5.1.1. The first animation corresponds to Example 13.3.11 with S = 2, a = b = 0, 
W(r) = V(r) = and the ip m (r,t) functions as in (|3.12p . namely (but with 3 different 2-vector 
parameters v m ) 

¥m (r, t) = A m + B m exp [-iv m (r - v m t)\ , 

with 



Ai = 1, A 2 = 5, A 3 = -5, 
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Figure 5.1.2. The absolute values of the ip m (x,t) functions and the absolute values of the ip m (x,t) 
solutions for < t < 20. 



Bi = 1, 

vi = (V5F,0), 



B 2 = 2, 

v 2 = (-2V5F,2>/37r), 



^3 = 3, 

v 3 = (-2 v / 7r,-2\/37r). 



Here the periods in time of the three ip m (r,t) functions are chosen to be the same, T = 1/8, 
and the animation is performed on the time interval < t < T, then closed in loop. 

Fig. 15.1.11 is the first frame. To see the whole animation, please click on the following 



(external) link: http:/ / www.emis.de/journals / SIGMA / 2006/Paper088 / Animation5. 1 .l.gif 



Example 5.1.2. The second animation corresponds to the Example 13.3.11 with S = 1, a = b 
0, W(r) = V(r) = and the ip m (x,t) functions as in (|3.14|) . namely 

\2" 



tp m (x, t) = A m + 



B, 



(t + t m + irj„ 



a/2 



exp 



i [x x m -\- lTl m V m ) 
4 (t + t m + ir] 

n 



with 





= 0, 


Bi 


= 0.5, 




= 30, 


Vl 


= 0.5, 


m 


= i, 


ti 


= 0.8 


A 2 


= -0.6, 


B 2 


= -0.9, 


x 2 


= 0, 


v 2 


= 1, 




= 2, 


t 2 


= 0.4 


A 3 


= 1.4, 


B 3 


= 0.8, 


X3 


= -20, 


V3 


= 0.5, 


m 


= 7, 


ts 


= 1. 



Here the animation is performed on the time interval < t < 20. 

Fig. 15.1.21 is the first frame. To see the whole animation, please click on the following 



(external) link: |http:/ /www.emis.de /journals /SIGMA / 2006 /Paper088 /Animation5. 1 . 2.gif 



5.2 Solutions of (TTT2I) 

In this subsection we present three numerical solutions of (jl.2p displayed as animations, two 
with S = 2 and the last with S = 1. 



Example 5.2.1. The first animation corresponds to the Example 14.3-H with S = 2, a = b = 0, 
M = and the <p m (x,t) functions as in (|4.4p . with K = 4 and a very particular choice of the 
functions ip m (x, y, t): 



ip m (x, y, t) = A m exp 



+ B m exp 



,(2) 



+ t 
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Figure 5.2.1. The <p m (r, t) functions and the absolute values of the ip m (r, t) solutions for < t < 20. 



where 



+ C m exp 



(!) + 
y-ym -t 



+ D rn exp 



y-ym +t 



+ E m , 



A! = -3, 


Si = 3, 


Ci = 0, 


£>i = 0, 


#1=0, 


A 2 = 0, 


B 2 = 0, 


C 2 = -6, 


D 2 = -3, 


E 2 = 10, 


A 3 = 6, 


£3 = 6, 


C 3 = 3, 


Z> 3 = 4.5, 


£ 3 = -12 


oi = 1.5, 


61 = 1.3, 


ci = 1.3, 


di = 1.6, 




a 2 = 1.8, 


62 = 1-2, 


c 2 = 2, 


d 2 = 1.6, 




a 3 = 1.4, 


63 = 1.5, 


c 3 = 1.4, 


d 3 = 1.2, 




*?> = -6, 


= 7, 


y? = 0, 


y ; 2) = 0, 




4 ] = o, 


4 2) = 0, 


4 1} = 2.5, 


y f = -1.5, 






4 2) = 6, 


y? = -5.5, 


yf ] = 4.5. 





Here the animation is performed on the time interval < t < 20. 

Fig. 15.2.11 is the first frame. To see the whole animation, please click on the following 



(external) link: http://www.emis.de/journals/SIGMA/2006/Paper088/Animation5. 2.1.gif 



Example 5.2.2. The second animation corresponds to the Example 14.3.21 with S = 2, a = b 
0, M = and the (p m (r, t) functions as in (j4.5f) : 



(x, y, t) = J ( y/(x- x m ) 2 + (y - y m ) 2 ) [A m cos(i) + B m sin(t)] + C m , 



where 



A 1 = 1, B 1= 0.1, Ci = 0, xi = 5, 



2/1 = 0, 
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Figure 5.2.2. The ip m (r, t) functions and the absolute values of the ip m (f, t) solutions for < t < 2ir. 

A 2 = -l, B 2 = l, C 2 = 3, x 2 = -5/2, y 2 = 5V3/2, 
A 3 = 0.1, B 3 = l, C 3 = -3, x 3 = -5/2, y 3 = -5Vs/2. 

The (fm (r,t) functions are periodic in time with period T = 2tt. The animation is performed 
on the time interval < t < T, then closed in loop. 

Fig. 15.2.21 is the first frame. To see the whole animation, please click on the following 
(external) link: |http:/ /www.emis.de/journals/SIGMA/2006/Paper088/Animation5.2.2.gifj 

Example 5.2.3. The third animation corresponds again to the Example 14.3. H with S = 1, 
a = b = 0, M = and the ip m (x, t) functions as in (|4.4p with K = 2, and with a very particular 
choice of the functions ip m (x, t): 

<p m (x, t) = A rn cos(x -t + B m ) + C m cos(x + t + D m ) + E m , 

where 

A 1 = 1, B x = 0, d = 3, D x = 0.1, E x = 0, 

A 2 = 2, B 2 = ir, C 2 = -l, D 2 =tt + 0A, E 2 = 5, 

A 3 = -l, S 3 = tt/2 + 0.1, C 3 = 2, D 3 = tt/2, E 3 = -5. 

The (p m (x,t) functions are periodic in time with period T = 2ir and periodic with respect to 
the x variable with period L = tt. The animation is performed on the time interval < t < T, 
then closed in loop. 

Fig. 15.2.31 is the first frame. To see the whole animation, please click on the following 



(external) link: http://www.emis.de/journals/SIGMA/2006/Paper088/Animation5. 2. 3.gif 



6 Outlook 

In future articles we plan to report additional investigations of these solvable (systems of) 
nonlinear evolution PDEs, and in particular to display visual animations of solutions of certain 
of these models in three-dimensional space. 
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Figure 5.2.3. The cp m (x, t) functions and the absolute values of the ip m (x, t) solutions for < t < 2tt. 
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